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Abstract. The behaviour of real zeroes of the Hurwitz zeta function 
C(s,a) =^(a + r)" s a>0 

r=0 

is investigated. It is shown that C( s i a ) nas no rea l zeroes (s = (J, a) in the region 
a > + log(— a) + 1 for large negative a. In the region < a < =^ the zeroes 
are asymptotically located at the lines a + 4a + 2m = with integer m. If N(p) is the 
number of real zeroes of £(— p, a) with given p then 

lim = J_. 

p— ioo p ne 

As a corollary we have a simple proof of Inkeri's result that the number of real roots 
of the classical Bernoulli polynomials B n (x) for large n is asymptotically equal to 22 _ 

1 Introduction. 

The classical Hurwitz zeta-function is defined for any positive real a as an analytic 
continuation of the series 

C( S ,a)=^> + rr. 

r=0 

When a = 1 it reduces to the Riemann zeta-function. 
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We should note that sometimes the definition of the Hurwitz zeta-function is re- 
stricted to < a < 1 (see e.g. Q). From our point of view this is not natural and 
we follow the definition of the Hurwitz zeta-function from J3j where all positive a are 
allowed (cf. also the original Hurwitz paper ^j). 

It is known (see e.g. volume 1, page 27) that in the special cases when s is 
a negative integer this function (as a function of the parameter a) reduces, up to a 
factor, to a Bernoulli polynomial: explicitly when s = -m, m = 0,1, 2, 3, .... 

st ™ ^ B m+ i(a) 



rn + 1 

The Bernoulli polynomials Bk(a) can be defined through the generating function: 



7 = £ 



Bk{a) u 



giving, for example: 



fe! 

fc=0 



1 „, s 2 .1 



B (a) = l Si(a)=a-- B 2 (a)=a-a+- B 3 (a) = a 6 - — + -, ... 

Bernoulli polynomials possess many interesting properties and arise in many areas 
of mathematics (see [[l], ^]). 

Inkeri || proved a remarkable fact that the number N(n) of real roots of the 
Bernoulli polynomials B n for large n asymptotically equals to — . More precise esti- 
mates for N(n) have been found by Delange || 

In this paper we investigate the behaviour of the real zeroes of the Hurwitz zeta 
function ((a, a) in the upper half-plane a > 0. As a corollary we have a simple proof 
of Inkeri's result. Our approach is different from ||, |^, Q and we believe is more 
elementary. It is based on the remarkable Hurwitz representation of the ((s,a) on the 
interval < a < 1 and Re(s) = a < 0: 

_ 2r(l - s) ^ sin(27rra + iyrs) 
C(S,a) ~ (2^)i- ~ ' 

r— 1 

where F is the Euler gamma- function (see e.g Q|, volume 1, page 26). 

Our main observation is that for a large negative a this formula gives a good 
approximation for ((a, a) on a much wider interval: < a < — As a result we 
prove that 



s^TT v = sin(27ra + -ira) + o(l) as a — > — oo provided < a < — 

Q(cr) 2 2ne 

where Q(s) = ^j^i-l ■ We show also that in the region a > — + log(— cr) + 1 
the Hurwitz zeta-function has no real zeroes for large negative a. 
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Figure 1 

Our results are illustrated in Figure 1 which shows the behaviour of the real zeroes 
of ((a, a). Since 7^ is a small number we have chosen different scales on the axes to 
make the picture more illuminating. Notice that £(s, 1) coincides with the Riemann 
zeta-function £(s) which only has real zeroes (for s < 0) at negative even numbers 
s = -2,-4,-6,.... Also, using the well-known identity £0,1/2) = (2 s - l)C(s), 
we see that the only real zeroes of C(s,a) on a = 1 are s = —2,-4,-6,... and on 
a = I they are s — 0, —2, —4, ... As we have mentioned above when s is a non-positive 
integer £(s, a) reduces to a Bernoulli polynomial. We have used this fact to compute 
numerically the corresponding values of £(s, a) and to draw the picture in the regions 
II and IV. When s = 1 the Hurwitz zeta-function has a simple pole with the residue 1 
and for a — Re(s) > 1 it is given by the convergent series with positive elements and 
therefore has no zeroes. 

2 Asymptotic behaviour of the Hurwitz zeta- 
function £(cr, a) for large negative a. 

The Hurwitz zeta-function (or generalised Riemann zeta-function) is defined as a series 
C(s,a) =^(a + r)~ s a>0 

r=0 
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in the complex domain Re(s) > 1 and can be analytically continued to a meromorphic 
function in the whole complex plane with the only pole at s — 1 (see (ffljfljPl). When 
a = 1 it reduces to the Riemann zeta-function 



C(*) = X>~ 
fc=i 



The Hurwitz zeta-function can be extended to the whole of the complex s-plane 
through the formula 

h ^ r(i- s ) f (0+) (_ z )-i e — , 

^^-^r 1 i- e - dz a>0 

</ oo 

in which the integral is taken over a curve starting at 'infinity' on the real axis, encircles 
the origin in a positive direction and returns to the starting point (see Q). By using 
an alternative integral formulation for ((s,a) it can be shown that a) is analytic 
everywhere except for the simple pole at s = 1. 

The Hurwitz zeta-function obviously satisfies the functional relation: 

n— i 

<(s,a) =((s,n + a)+J2(r + a)- 3 n = l,2,... (1) 

r=0 

Since each term in this relation is analytic we can assume this relation is true for the 
whole of the complex s plane, except for s = 1. 

In this paper we restrict ourselves to the case when s is real: s = a £ R. When a 
is negative Hurwitz has found the following Fourier representation for £(<j, a) on the 
interval < a < 1: 



(2tt) 1 - ct ^ r 



r=l 



From this formula we see that 



C(<r, a) 1 . 

= sm(27ra + — ira) + o(l) when o —> —oo, 

Q(<t) 2 

where Q(u) — 2 ^\i°} ■ Our first theorem proves that this is actually true on a much 
larger interval. 

As part of the theorem proofs we will use the following inequality for the function 
S(p,n) = V 1 +2 p + ... +n p : 

S(p, n) < n p 



Indeed, 

S(p, n) = n p (l + (1 - -Y + (1 - -Y + . . . + (1 - — -Y ) 

\ n n n J 

But since 1 — x < e~ x for x < 1 we have 

1 ^ ^ -l/n 1 ^ _2/n -| n— 1 ^ -(n-l)/n 

1 <e' 1 < e ' ... 1 < e K " 

n n n 
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and therefore 

S(p,n) <n p (l + e- p/n + e- 2p/r ' 



+ ■ 



- (n—l)p/n 



)=»' T 



l-e- p 

-p/n 



We shall also need an estimate for the sum of the series: y^°° — i- . Now 



11 1 1 fl .1,2 

2 1 +p 3 1+ p 4!+p ' ' " — 2 p 



V2 3 V 3 ; 4 V 4 ; / 



But < i if p > 4 and n > 3 so we easily deduce 



oo 

E < («2) " |)2" P P > 4 



Finally, we will also make particular use of Stirling's inequality for the gamma function: 
(2np)?p p e~ p < r(l +p) < (2-Kp)ip p e~ p e T ^ p » 1 

Theorem 1. Let a — —p,p > and < a < ap for some positive a then the Hurwitz 
zeta-function satisfies the inequality 



Q(-p) y 2 L} 



<C lp - 1/2 (2nea) p + C 2 2- 



(3) 



where C\,Ci are constants, which do not depend on p. In particular, on the interval 
< a < 2^P we have the asymptotic behaviour 



^ -5**) + °(D 



when p — > oo. 



Proof 



Let us represent aasn + 6, < & < 1 and with n integer. It follows from the 
functional relation (|l|) that 



C(a,a) C(v,b) 



Q(a) Q(a) 



C(a,n + &) C{a,b) 



Q(a) Q(o-) 



n-l 

^E^r 



Since < b < 1 we obviously have 



^(r + &)- CT <S(p, 



< p = -cr 



and, as we obtained above, 



S{p, n) < n p 



1 - e- p 
1 - e -p/" 



Also, from the Stirling formula for the F-function we have the following asymptotically 
exact inequality T(l + p) > (2np) 1 ^ 2 p p e~ p and therefore 

1 (2tt) 1+p r 27re 1D tt 

< { — Y 



Q(-p) 2r(l+p) P 1 */EFp~ 
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Thus for a large p 

C(-P,a) C(-P,b) 



Q(-p) Q(-p) 



I 2Tven 



P 

2-nen 



1 - e" 



2np \ 1 — e 



/Sip (1 



-p/n 



if - < a 



(4) 



which is true since ~ < ~ < a by assumption. 
However, from Hurwitz' formula (^) it follows that 

Cf— P, 6) sin(27rrfc — i7rp) 1 sin(47rfc — i7rp) sin(67rfo — ^?rp) 
— - - — 1+p 2 = sm(27r6--7rp)H ^J—^ 



Q(-p) 

Therefore 

a-pM 



2 i+p 



-+- 



3 i+ P 



+. 



q(- P ) ^ h -- 2 - p) 



a-p,b) . .„ 1 , 



r = 2 



(5) 

if p > 4. The estimates (^) and (|^) imply the theorem. 

By a slight modification of the previous arguments we can prove the following result. 
Theorem 2 In the region 

a> £-e + 4k l ° gP+1 
f(— p,a) is negative if p is sufficiently large. 
Proof From the same functional relation ([!]) 

C(-p,q) c C(-p,ft) 
0(-p) Q(-p) Q(-p) ' 

Now assuming that p is sufficiently large we can use Stirling's inequality T(l + p) < 
v /27rp(f) P e T ^, so 

However, we know from the Hurwitz formula, that when p — > 00 
^^= s in(27r6-i7rp) + (l) 

Therefore if ^V^L > 1 and p large enough then £(— p, a) < 0. But as we have shown 



Q(-p) 



(o-lf / 27re(o-l) y /V 



0(-p) > I 



2p 



-e 12 f 



which is greater than f if 



p /2p\2p _1_ p { 



27re \ 7r 



(log 2p-log tt) , 1 

2p ~ l ~12? 3 '' 



2vre 



Now, using the inequality e x < 1 + x + x 2 for sufficiently small x we find that to 

guarantee that > 1 it is enough to demand that 

Q(-p) to 



This implies the theorem. 



P 1 , 

27re 47re 



G 



3 The Real Zeroes of the Hurwitz Zeta-function 
and Bernoulli Polynomials. 

Let us now fix a — — p and consider C(~Pj a ) as a function of a. It follows from 
theorem 2 that the zeroes of this function for large p are located in the interval < 
a ^ 2yre + ike 1°8 P + 1- F° r given p let N(p) be the number of real zeroes of C(~P> a )j 
and A(p) be the largest of these zeroes. 
Theorem 3 For p sufficiently large 

1 <Mp)<& + -±-Jo gP +l (6) 



27re 2 27re 47re 

— - - 1< N(p) < — - + ilogp + 27re + 2 (7) 
ne Tie 2 

The zeroes of f (— p, a) on f/ie interval < a < 2=- are simple and close to the 

half-integer lattice: a = 2 + i i 6 Z. 

Proof Let us introduce the function Z„(a) — < *k7 p ' a ? . From theorem 1 it follows that 

py ' Q(-p) 

Z p (a) = sin(27m — ^ 7r P) + 

on the interval I p : < a < when p — > oo. Actually this is true also for the A; th 
derivative of Z p (a) but on a smaller interval I p -k- Indeed, from the definition of the 
Hurwitz zeta-function we see that 

^C(s,a) = (s)C(s+l,a). 

From the property of the F-function F(p + 1) = pF(p) it follows that 

«-p>-^-£«-p + «- 

Fhus the derivative of Z p (a) is equal to 

1 1 

Z p (a) — 2n Zp-i(a) — 2n sin(2-7ra — ^(p — 1)) + o(l) = 27r cos(27ra — T^V) + "(1) 

on the interval J p _i. Similarly we have for the fe th derivative of Z p (a) 
Z p k) (a) =sin (fe) (27ra- ~7rp) + o(l) 

on the interval I p -k- 

In particular, on the interval I p -i the function Z p (a) (and its derivative) tend 
to sin(27ra — ^np) (and its derivative) when p — > oo which ensures that for large 
p all the roots of £(— p, a) on this interval are simple and located near the points 
ffl=f + 5 I E Z. This implies the last statement of the theorem and the lower 
estimates of (|^) and (Q). 

The upper estimates for A(p) follows directly from theorem 2. To prove the upper 
estimates for N(p) we need the following simple lemma. 

Lemma If a function f(x) (with a continuous n th derivative) on some interval (a, 6) 
has the property that the sign of the n th derivative is constant throughout the interval 
then / has no more than n roots on this interval. 
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Now we apply this lemma to the function Z v (a) on the interval J p : (5—^-, + 
t^t logp + 1) to estimate the number of roots there. The idea of this calculation is 
clear from Figure 2 (in which k = 5=-). 




Figure 2 

Using the fact that 

4 n) (a) = (2it) n Z p - n (a) 
we differentiate Z p (a) many times until we have a negative function and then apply 
the lemma. As one can see from Figure 2 if n > [y] then Zp (a) will be negative in 
the interval J p and, as such, cannot have more than n simple roots in this interval. 
Now y is the solution to the equation 

« Up - v) + g 1o s(p - y)) + 1 = k (p - !) 

or 

-y+-log{p-y) + (K- 1 + 1) = 
We claim that the solution to this equation for large p satisfies the inequality 

y < - logp + 2ne + 1 

Indeed the function F(y) = — y + | log(p — y) + {2ixe + 1) is monotonically decreasing 
and 

F{- logp + 27re + 1) = --logp- 2ixe - 1 + ~log(p- i logp - 27re - 1) + (27re + 1) 

1 i ci 1 i 27re+l. ^ n 

2 2p p 

for large p. Thus according to the lemma Z v (a) has no more than | logp + 27re + 1 
roots on the interval (%r—, tt — h -r— logp+ 1). Since on the interval (0, %r—\ we have 
no more than — — h 1 zeroes this proves theorem 3. 

As we have already mentioned, when p = — a — m, m € Z + the Hurwitz-zeta 
function reduces to certain polynomials related in a simple way to the Bernoulli poly- 
nomials: 

a/ \ B m+ i(a) . . 
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Theorem 3 applied to these special values of p gives some estimates on the real pos- 
itive roots of the Bernoulli polynomials but because of their well-known symmetry 
properties: 

flm(l-a) = (-l) m B m (a) (9) 

we can immediately extend this result for all real roots of B m (a). In particular if 
N(m) is the number of all real roots of B m (a) and A(m) is the largest of these roots 
then from Theorem 3 it follows that for a large m 

m fl M .,n m 1, i, 1 » 

h - < Mm) < 1 logm + 1 (10) 

2ire \ixe 2J K ' 2ne Ane 6 ^ 2k e J K ' 

— -( —+2] <N(m)< — +logm+ (iire + 3-— ^ (11) 
7re \ ire J ire \ ne I 

Corollary (K. Inkeri). 

N(m) 2 , A(ml 1 

lim = — , km = - — . 

m— >oo to 7re m^oo m 2ne 

Remark. H. Delange in || ^ has found sharper estimates for A(m) and N(m). 
In particular he showed that the additional logarithmic terms exist in both upper 
and lower bounds. This should be true also for the real zeroes of the Hurwitz 
zeta-function but it does not follow from our elementary arguments. 
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